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We report the results of a search for the anomalous production of a massive particle decaying to 
four electrons via two bosons in 1.1 fb~^ of pp collisions at ^/s =1.96 TeV collected by the CDF II 
detector at Fermilab. We employ optimized electron identification criteria to maximize acceptance 
and efficiency. We estimate the backgrounds in the invariant mass range 500 - 1000 GeV/c^ to be 
0.028 ±0.009 (stat) ± 0.011 (syst) events. We observe zero events in this search region. Assuming a 
Randall-Sundrum graviton production model, we set 95% CL limits on cr x BF(G^ Z°Z°) < 4-6 
pb, depending on the graviton mass. 



PACS numbers; 12.60.Cn,13.85.Rm,14.70.Hp 
I. INTRODUCTION 



We present a search for new heavy particles 
"G" in the decay mode G ^ Z'^ Z'^ eeee in pp 
colhsions at -^/s = 1.96 TeV performed with the 
CDF II detector at the Fermilab Tevatron. Pre- 
vious analyses of double gauge boson production 
at the Tevatron have focused on standard model 
production [H [2 [3l H] . Here, we present for the 
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first time a search for massive particles G which 
decay to Z'^Z'^ which could indicate physics be- 
yond the standard model. 

The goal of this search is to be sensitive to 
production of any massive particle which could 
decay to Z°Z°, that is, to avoid focusing on any 
one specific model; however, for the purpose of 
quantifying acceptance for this search, we con- 
sider the virtual production of gravitons in a 
simple Randall-Sundrum (RSI) scenario [SJ 
In this model, the geometry consists of two 3- 
branes which confine the standard model sector 
separated from each other by a single extra di- 
mension. One can look for evidence of the ex- 
tra dimension at particle colliders in the form of 
a Kaluza-Klein tower of discrete, massive gravi- 
tons. In the RSI scenario, the gravitons predom- 
inantly decay to jets, and the remaining modes 
are W+W- (~ 10%), Z°Z° 5%), 77 5%), 
and II 2% per lepton)[7, . Searches for the de- 
cays of such particles to photons and electrons 
have been performed [HI M US]- If the couplings 
to leptons and photons are suppressed relative to 
the couplings to gauge bosons [U^ , such a particle 
might escape detection in these searches. Here, 
we have searched for massive particles in their 
decays to Z° bosons. 

In the leptonic final states of Z° decay, the 
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expected signal from the model described above 
is small, as are the backgrounds. In order to 
maximize acceptance and efficiency for the four- 
electron signature, we use optimized calorimet- 
ric electron identification criteria, select electron 
candidates identified as isolated tracks where 
there is no calorimeter coverage, and use low elec- 
tron energy thresholds. 

The organization of this article is as follows: 
first, we describe the components of the CDF II 
detector relevant to this search and summarize 
the data sample and event selection criteria. 
Then we describe the background estimation, re- 
port the results of the search, and interpret the 
results in the context of the lightest massive RSI 
graviton. 



II. THE CDF II DETECTOR 

This analysis uses 1.1 fb""'^ of pp collisions col- 
lected by the CDF II detector, a general purpose 
magnetic spectrometer. We briefiy describe the 
components of the detector relevant to this search 
here. A complete description can be found else- 
where [12]. 

A combination of tracking systems reconstruct 
the trajectories and measure momenta of charged 
particles in a 1.4T solenoidal magnetic field. Tra- 
jectories of charged particles are reconstructed 
using an eight-layer silicon microstrip vertex 
tracker |T3] at radii 1.3 < r < 29 cm |T3], and 
a 96-layer open-cell drift chamber (COT) pro- 
viding eight superlayers of alternating axial and 
stereo position measurements [TS] . The COT al- 
lows track reconstruction at large radii 43 < r < 
132 cm in the region I77I < 1.6, and provides full 
geometric coverage for I77I < 1.0. 

Outside the tracking volume, segmented elec- 
tromagnetic (EM) lead-scintillator and hadronic 
(HAD) iron- scintillator sampling calorimeters 
measure particle energies [16j . In the cen- 
tral region (jTyj < 1.1), the calorimeters are 
arranged in a projective-tower cylindrical ge- 
ometry, divided azimuthally into 15° wedges 
which measure EM energies with a resolution 
of HE)/E]'' = {13.5%)yET + (2%)2. In 
the region covered by the forward calorimeter 
(1.1 < \ri\ < 3.6), the calorimeters are ar- 
ranged in an azimuthally-symmetric disk geome- 
try and measure EM energies with a resolution of 
HE)/E]'' = (16.0%)V£;+ (1%)2. Wire cham- 
bers (scintillator strips) embedded in the cen- 
tral (forward) EM calorimeters at the electro- 



magnetic shower maximum QXq) provide po- 
sition and lateral shower development measure- 
ments used to identify electrons by their char- 
acteristic energy-deposition distribution. The 
beam luminosity is determined by measuring the 
inelastic pp collision rate with gas Cherenkov 
detectors [T7], located in the region 3.7 < < 
4.7. 



III. EVENT SELECTION 

Events are selected for collection by a three- 
level trigger system. We search in data collected 
by triggering on a central high-momentum elec- 
tron. Each of two trigger paths used requires an 
energy cluster in the central calorimeter and a 
track which projects to the energy cluster. The 
primary trigger path requires a clustered trans- 
verse energy Et > 18 GeV, transverse momen- 
tum of the associated track px > 9 GeV/c, the 
ratio of energy measured in the hadronic to elec- 
tromagnetic calorimeters, Ehad/ Eem < 0.125 
and a lateral shower profile consistent with an 
electron. The second trigger path requires a clus- 
ter with i?T > 70 GeV and an associated track 
with pt > 15 GeV/c. 

We select events containing one "seed" electron 
which satisfies trigger requirements and those 
listed in Table |T] and three other electrons which 
satisfy either the central, forward, or track re- 
quirements in Table |l] To maximize acceptance 
and efficiency for events containing four elec- 
trons, we select three other electron candidates 
using optimized identification and kinematic cri- 
teria in the central or forward calorimeters, and 
from isolated tracks pointing to uninstrumentcd 
regions of the calorimeters. 

Electron candidates are formed in the central 
and forward calorimeters from isolated energy de- 
posits with Et > 5 GeV. An electron is con- 
sidered to be isolated in the calorimeter if the 
sum of the transverse energy detected within 
a cone AR = y/iAr])"^ + (A^)^ < 0.4, minus 
the electron Et, is less than 20% of the elec- 
tron Et. For clusters in the central calorime- 
ter, where tracking efficiency is optimal, we re- 
quire that a track reconstructed in the COT 
project to the cluster. Tracks must include mea- 
surements in at least three axial and two stereo 
superlayers of the COT, and the track coordi- 
nate along the beam direction, zq, must also 
lie within the nominal extent of the interac- 
tion region, \zq\ < 60 cm. To reduce back- 
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grounds from hadrons misidentified as electrons, 
central candidates must also satisfy an energy- 
dependent requirement Ehad/Eem < 0.055 + 
0.00045 GeV~^ x Eem- Forward candidates must 
have Ehad/Eem < 0.05 and he within 1.1 < 
hi < 2.5. 

The "seed" electron candidate must satisfy the 
above requirements to be reconstructed in the 
central calorimeter, and satisfy additional selec- 
tion criteria imposed by the trigger. Specifically, 
the "seed" electron must also have Et > 20 GeV, 
satisfy the same lateral shower profile require- 
ment as the triggering one, and have an associ- 
ated track with px > 10 GeV/c. 

Some calorimeter acceptance for electrons is 
lost in 24 one-degree gaps in (j) between the cen- 
tral calorimeter projective wedges, a region at 
0.7 < rj < 1.0 and 75° < cj) < 90° which 
accommodates cables and cryogenic utilities for 
the solenoid, a gap between central calorimeter 
arches at = 0, and the gap at the junction 
between the central and forward calorimeters at 
1.0 < \r]\ < 1.2. Together, these regions add 
up to approximately 17% of the solid angle for 
I77I < 1.2, which for a four-lepton final state rep- 
resents an acceptance loss of approximately half. 

We recover acceptance by forming electron 
candidates from isolated tracks which project to 
the gaps between instrumented regions of the 
calorimeter. A track is defined to be isolated in 
the tracking chamber if the transverse momen- 
tum of the track is more than 90% of the total 
transverse momentum of all tracks within a cone 
AR < 0.4 around the candidate track. We re- 
quire track electron candidates be consistent with 
originating from prompt decays by requiring that 
they pass within 0.2mm (2mm) of the axial beam 
position for tracks with (without) position mea- 
surements in the silicon vertex tracker. 

To reconstruct Z° candidates, we form all 
unique combinations of pairs of electrons in the 
event. To avoid rejecting events where the charge 
of one electron is misidentified, we do not impose 
an opposite charge requirement on the pair. We 
ensure that all electron candidates are isolated 
from each other by requiring a separation of 0.2 
in AR between any two electron candidates in 
the combination. If both candidates in a pair 
have associated tracks, we ensure they are con- 
sistent with originating from the same parent by 
requiring their zq measurements to lie within 5 
cm of each other. The invariant mass distribu- 
tions for events containing Z*^ candidates formed 
from a "seed" electron candidate together with 




50 100 150 

(GeV/c") 

FIG. 1: Distribution of mee for Z'^ candidates formed 
from a "seed" electron candidate together with a sec- 
ond electron candidate (a), and the subset of can- 
didates formed from a "seed" electron candidate and 
an isolated track (b). 



just one other electron candidate and the sub- 
set where an isolated track is used as the second 
electron candidate are shown in Fig. [T] 

To reconstruct Z'^ pairs, we form all unique 
combinations of all Z° candidates containing a 
"seed" electron with all other candidates in 
the event and again require a separation of 0.2 
in Ai? between any two electrons in the four- 
electron combination. 

The variable 




quantifies consistency between a given combi- 
nation and a Z'^Z''^ —> eeee final state, where 
— 91.19 GeV/c^ is the nominal Z° mass [TS], 
rrii is the measured invariant mass of each candi- 
date Z'^ in the combination computed from the 
electron candidates' four-momenta, and ai is the 
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TABLE I: Electron Selection Criteria 



Criteria 



Central (Seed) 



Forward Track 



Et (GeV) 
1 20 1 (cm) 
Ehad/Eem < 0.055 

Isolation 
PT (GeV/c) 



> 5(20) 

< 60 
0.00045 GeV" 

< 0.2 
(> 10) 



X E 



> 5 

< 60 

< 0.05 

<0.2 >0.9 
> 10 



uncertainty on the mass of each Z'^ candidate 
consisting of a contribution from the intrinsic 
width of the Z'^ boson and a contribution propa- 
gated from the individual electron energy or mo- 
mentum measurements (typically ^ 3.5 GeV/c^). 
In each event, we retain the one Z'^Z'^ combina- 
tion with the lowest x^. 

We fix the final event selection criteria before 
examining the event yield in the signal region. 
We define the signal region as the events con- 
taining a four electron combination with < 50 
and nieeee > 500 GeV/c^. Events with rrieeee < 
400 GeV/c^ are used as a low-mass control re- 
gion to estimate backgrounds, as described in 
Section |IV[ We find 12 events containing four 
electron candidates in this control region. 

Although we do not focus on one spe- 
cific model, for the purpose of quantifying 
acceptance for this signature, we consider a 
graviton-production scenario implemented in the 
HERWIG [20] Monte Carlo generator which is 
treated in a model-independent way, assuming 
only that there is a universal coupling of the 
graviton to standard- model particles. For com- 
parison, we interpret the couplings in the con- 
text of the RSI model. We determine geomet- 
ric and kinematic acceptance and reconstruc- 
tion efficiency for this model using Monte Carlo 
calculations followed by a GEANT-based simula- 
tion of the CDF II detector [T^. We consider 
graviton production followed by decay to two 
bosons followed by decay into four electrons. We 
use a leading-order calculation implemented in 
HERWIG to estimate acceptance times efficiency 
for the model. For a RSI graviton with mass 
Mg = 500 GeV/c^ and ratio of warp factor to 
Planck mass, k/Mpi = 0.1, Fig.[2]shows the dis- 
tribution of reconstructed ^^id 

TTZgeeei the in- 
variant mass of the Z'^Z'^ combination with the 
lowest computed from the four-momenta of 
the two Z° candidates. As expected for events 
containing two real bosons, the x^ distribu- 
tion peaks near zero, and the total invariant mass 
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FIG. 2: Distribution of for simulated Randall- 
Sundrum signal scenario {ma = 500 GeV/c^) (top). 
Four-electron invariant mass distribution for events 
satisfying < 50 (gray) and for events which fail 
this requirement (black) (bottom). 



of selected combinations is centered on the gener- 
ated graviton mass (500 GeV/c^). Events which 
contain mis-measured electrons contribute to the 
population with large x^ values. The width of 
the rrieeee distribution, 15 GeV/c^, is dom- 
inated by the detector resolutions of the con- 
stituent electron candidates. We find the geo- 
metric and kinematic acceptance times efficiency 
for this model to be 65%. Of the events recon- 
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FIG. 3: Acceptance for Randall-Sundrum graviton 
decaying to Z'^Z'^ versus its mass. 



structed, 93% yield a four electron combination 
with < 50. 

The total acceptance versus rrieeee is shown in 
Fig. |3] At very high graviton mass, the mo- 
mentum of the daughter Z^s becomes significant, 
which can cause the electrons to have a small 
opening angle and fail the isolation requirement. 



IV. BACKGROUND ESTIMATION 

In studies using Monte Carlo simulation to es- 
timate the main sources of backgrounds, we find 
that the dominant background consists of events 
in which one or more hadrons satisfy the electron 
ID requirements in the four electron combination. 

We use control samples in the data to obtain 
the shape and normalization of this background 
in the signal region. Background-dominated 
(hadron-enriched) control samples are selected 
from the data. We form hadron candidates, h, 
from calorimeter clusters in a manner identical 
to the central and forward electron candidates, 
with two exceptions. The hadron candidate must 
fail the relevant Ehad/ Eem criterion, and to 
increase the size of the control samples, we do 
not impose any isolation requirements. In Fig. [4] 
we show the invariant mass of all pairings of one 
seed electron candidate with one hadron candi- 
date. The absence of a significant peak at the Z'^ 
mass indicates that contamination from electrons 
in the hadron candidates is small. ^ 



^ The presence of a small amount of electron contami- 



6000^ 



> 
a 

•c 




50 100 150 

m^, (GeV/c') 

FIG. 4: Invariant mass distribution of one cen- 
tral electron satisfying trigger requirements and one 
hadronic candidate in data. 



We obtain five control samples, namely 
the four-electron sample which has nieeee < 
400GeV/c^ introduced above, and additional 
control samples having one, two, three, or four 
hadron candidates by repeating the Z^Z^ selec- 
tion procedure, forming combinations using one 
or more hadron candidates with electron candi- 
dates and retaining the minimum combination 
for each event. The distributions of the minimum 
versus rrieeee for samples with different num- 
bers of hadron candidates are shown in Fig. [5] 
For reconstructed masses smaller than twice the 
mass, there is a correlation between ^^^d 
mass caused by a kinematic threshold effect. At 
higher masses, the two variables are much less 
correlated. 

The single probability density function, 

fiX^, rrieeee) = Cm^eeeC^'^, (2) 

where C is a normalization constant, provides an 
empirical description of the x^ vs. rueeee dis- 
tributions for each of the four hadron-enriched 
control samples. We obtain the parameters 7 — 
-4.57 ± 0.06 and r = -0.00319 ± 0.00007 from 
a two-dimensional unbinned maximum likelihood 
fit to the low-mass four-electron control region 
and the hadron-enriched control samples simul- 
taneously, using events with invariant mass above 



nation in the hadron candidate sample has a negligible 
effect on the estimate of the background at high mass, 
and is included in the systematic uncertainty we assign 
to the background estimation method. 
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FIG. 5: Distribution of vs nieeee for control samples containing one, two, three, or four hadron candidates 
with the number of events in the plot increasing with the number of hadron candidates used in the combination. 



185 GeV/c^ 2 X rn^o.) The control samples 
containing mostly hadron candidates dominate 
the fit. Fig.|6]shows the projections of the fit re- 
sult in the invariant mass dimension along with 
the data for each hadron-enriched sample. The 
low-mass control region {nieeee < 400 GeV/c^) 
contains five events with meeee > 185 GeV/c^ 
and serves to normalize the prediction of back- 
ground in the high-mass search region. We inte- 
grate the fit result above 500 GeV/c^ and < 50 
to extract an estimate of the total background 
for the high-mass region. Using this method, we 
estimate 0.020 ±0.009 (stat) ± 0.007 (syst) back- 
ground events from hadrons misidentified as elec- 
trons in the search region. The systematic uncer- 
tainty on the background estimate is obtained 
by varying the functional form of the probability 
density function fitted to the rUeeee spectrum. 

We have performed several cross checks to en- 
sure that the fit provides a reasonable estimate of 
the background rate at high mass. In particular, 
we have performed the fit allowing the power-law 
parameter to vary independently for each of the 



TABLE 11: Result of fit with 7 floating independently 
for each control sample. 



Sample 


Events 


7 


eeee 


5 


-5.90 ±2.14 


eeeh 


52 


-4.85 ±0.55 


eehh 


323 


-4.28 ±0.19 


ehhh 


1208 


-4.43 ±0.10 


hhhh 


1927 


-4.71 ±0.09 



categories. The parameters resulting from this 
fit along with the number of events observed in 
each sample are shown in Table |ll] The fitted 
values for 7 are consistent within errors across 
categories and with the nominal result. We have 
verified that the background shape in meeee is in- 
dependent of in subsets of data in bins of 
and have checked that the projected fit result is 
consistent with data. 

Standard model production of Z^Z'^ [21 is the 
only background to this search with two real 
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FIG. 6: Projections of fit to invariant mass in control samples of varying number of electron and hadronic 
candidates with same ordering as in Fig. [5] Data are shown with the fit projection overlaid. 



bosons and possibly four electrons in the final 
state. While we use data to estimate the total 
background from misidentified electrons in the 
search region, we have studied the background 
from this source at high mass with simulated 
events. We determine geometric and kinematic 
acceptance using Monte Carlo events generated 
by PYTHIA , followed by a GEANT-based sim- 
ulation of the CDF II detector. The expected 
number of events from this background compo- 
nent is determined as the product of the cross 
section, the luminosity of the sample, and the 
acceptance of the detector. We estimate a total 
of 0.54 ± 0.04 events in the four-electron sample. 
In the invariant mass region above 500 GeV/c^, 
we expect 0.008 ± 0.006 events. We estimate the 
total background from production of standard 
model Z'^Z'^ events and events in which hadrons 
are misidentified as electrons is 0.028 ± 0.009 
(stat) ± 0.011 (syst) events. 



OTHER SYSTEMATIC 
UNCERTAINTIES 



When setting the cross-section limit, we have 
considered other systematic uncertainties from 
several sources. The dominant source of these 
uncertainties is from the measured luminosity 
(5.9%) p3]. Other sources include parton dis- 
tribution function uncertainties (0.4%), signal 
Monte Carlo statistics (1.3%), initial state radia- 
tion(1.0%), and the difference between electron 
identification efficiency in data and simulation 
(1.0% per electron). The total systematic un- 
certainty from all these sources is 7.3%. 



VI. RESULTS 

The distribution of data events surviving all 
requirements is shown in Fig. [7] We observe no 
events in the high-mass signal region. There is 
one event in the low-mass region with very small 
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FIG. 7: Distribution of x'^ vs meeee for four-electron 
candidates in data. The signal region is < 50 and 
meeee > 500 GeV/c^, the 
shown by the dashed line. 



background from standard model processes in- 
volving four electrons with simulations. Us- 
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boundaries of which are piQ. 8: Limits on cr x BF(G ^ Z° Z^) versus mass. 



consistent -with standard model Z'^Z^ pro- 
duction, "while we expect 0.54 standard model 
Z^Z" events over the entire mass range. In 
the observed event, the total invariant mass is 
190 GeV/c^, and the two Z^ candidates in the 
lowest combination have measured masses of 
91 and 92 GeV/c^. 

We have set hmits on a(j)j) ^ G) x BF(G 
Z'^Z'^) in the context of a RSI graviton scenario. 
We use a Bayesian binned maximum likelihood 
method to extract 95% confidence level limits 
in 100 GeV/c^ wide windows centered on each 
mass. We incorporate the effects of uncertainty 
on the background and signal acceptance with a 
flat prior for the signal rate and Gaussian pri- 
ors for the acceptance and expected background. 
The limits including systematic uncertainties on 
the acceptance on ctxBF(G Z^Z^) range from 
4— 6 pb, depending on the graviton mass, and are 
shown in Fig. [8] along with the prediction from 
the RSI model for k/Mpi = 0.1. In the future, 
the sensitivity of this search will improve with 
more data as well as additional acceptance from 
other Z° decays. 

VII. CONCLUSIONS 

We have searched for production of particles 
decaying to a pair of bosons. We have 
estimated backgrounds from misidentified elec- 
trons using a data-based technique, and the 



ing an optimized electron selection, we expect 
0.028 ± 0.009 (stat) ±0.011 (syst) total back- 
ground events with < 50 above 500 GeV/c^ 
in 1.10 fb~^, and observe no events. In the ab- 
sence of evidence for a signal, we have set limits 
on G(1.96GeV)) x BF(G 



Acknowledgments 

We thank the Fermilab staff and the technical 
staffs of the participating institutions for their 
vital contributions. This work was supported by 
the U.S. Department of Energy and National Sci- 
ence Foundation; the Italian Istituto Nazionale di 
Fisica Nucleare; the Ministry of Education, Cul- 
ture, Sports, Science and Technology of Japan; 
the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council 
of the Republic of China; the Swiss National Sci- 
ence Foundation; the A. P. Sloan Foundation; the 
Bundesministerium fiir Bildung und Forschung, 
Germany; the Korean Science and Engineering 
Foundation and the Korean Research Founda- 
tion; the Science and Technology Facilities Coun- 
cil and the Royal Society, UK; the Institut Na- 
tional de Physique Nucleaire et Physique des Par- 
ticules/CNRS; the Russian Foundation for Basic 
Research; the Comision Interministerial de Cien- 
cia y Tecnologi'a, Spain; the European Commu- 
nity's Human Potential Programme; the Slovak 
R&D Agency; and the Academy of Finland. 



12 



Phys. 



Phys. 



[11 A. Abulencia et al. (CDF Collaboration 

Rev. Lett. 98, 161801 (2007). 
[21 V. M. Abazov et al. (DO Collaboration; 

Rev. Lett. 95, 141802 (2005). 
[3] A. A. Affolder et al. (CDF Collaboration), Phys. 

Rev. Lett. 88, 071806 (2002). 
[4] D. E. Acosta et al. [CDF Coll aboration], 
Phys. Rev. D 71, 091105 (2005) [arXivihep- 



[9] 
[10] 

[11 

[12] 
[13] 

[14] 



[15] 

[16] 

[17] 

[18] 
[19] 

[20] 
[21] 



ex/050102H. 
L. Randall and R. Sundrum, Phys. Rev. Lett. 
83, 3370 (1999). 

L. Randall and R. Sundrum, Phys. Rev. Lett. 
83, 4690 (1999). 

B. C. AUanach, K. Odagiri, M. J. Palmer, 
M. A. Parker, A. Sabetfakhri and B. R. Web- 
ber, JHEP 0212, 039 (2002) 
T. Aaltonen et al. [CDF Collaboration], Phys. 
Rev. Lett. 99, 171801 (2007) 
T. Aaltonen et al. (CDF Collaboration), 



|arXiv:h ep-ex/07072521. 

(DO Collaboration) 



Phys. 



V. M. Abazov et al. 
Rev. Lett. 95, 091801 (2005). 
A. L. Fitzpatrick, J. Kaplan, L. Randall, and 
L. T. Wang, J. High Energy Phys. 0709, 013 
(2007). 

D. Acosta et al. (CDF Collaboration), Phys. 
Rev. D 71, 032001 (2005). 

A. Sill et al, Nucl. Instrum. Methods A 447, 
1 (2000). A. A. Affolder et al, Nucl. Instrum. 
Methods A 453, 84 (2000). 

CDF uses a cylindrical coordinate system in 
which 9 {(f>) is the polar (azimuthal) angle, r is 
the radius from the nominal beam axis, and +z 
points along the proton beam direction and is 
zero at the center of the detector. The pseudo- 
rapidity is defined as rj = — lntan(&/2). Energy 
(momentum) transverse to the beam is defined 
Et = Esind (px = psind), where E is energy 
measured by the calorimeter and p is momentum 
measured by the spectrometer. 
A. A. Affolder et al., Nucl. Instrum. Methods A 
526, 249 (2004). 

L. Balka et al., Nucl. Instrum. Methods A 267, 
272 (1988). 

D. Acosta et al, Nucl. Instrum. Methods A494, 
57 (2002). 

W.-M. Yao et al, J. Phys. G 33, 1 (2006). 
R. Brun, F. Bruyant, M. Maire, A. C. McPher- 
son and P. Zanarini, "GEANT3," CERN Report 
No. CERN-DD/EE/84-1 (unpublished). 
G. Corcella et al., "HERWIG 6: An event gen- 
erator for hadron emission reactions with in- 
terfering gluons (including supersymmetric pro- 
cesses)," J. High Energy Phys 0101, 010 (2001). 
J. M. Campbell and R. K. Elhs, Phys. Rev. D 
60, 113006 (1999). 



[22] T. Sjostrand, S. Mrenna, and P. Skands, 
"PYTHIA 6.4 physics and manual," J. High En- 
ergy Phys. 0605, 026 (2006). 

[23] S. Klimenko, J. Konigsberg, and T.M. Liss, 
FERMILAB-FN-0741 (2003) (unpublished). 



